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1. Introduction 

Lattice QCD calculations of the properties of hadronic interactions such as elastic scattering 
phases shifts and scattering lengths have recently started to develop. Precision results have been 
obtained in the meson-meson sector for certain processes such as pion-pion, kaon-kaon and pion- 
kaon scattering and preliminary results for baryon-baryon scattering lengths have been presented. 
A recent review of these calculations can be found in[l]. In this work we study scattering processes 
where one hadron contain a charm quark. Firstly, we study the scattering processes of chamonia (t]c 
and with light hadrons (;r, p and N). As it has been pointed out in the literature [2, 3, 4], such 
interaction has a direct relation to possible charmonium-nucleus bound states with binding energy 
of a few MeV. Unlike the traditional nuclear force that binds nucleons, in this case, there are no 
quark exchange diagrams, and only gluons are responsible for the binding. In other words, the 
charmonium nucleon force is purely a gluonic van der Waals force. The charmonium interactions 
with light hadrons has been also studied in quenched lattice QCD [5]. Secondly, we study the 
scattering processes of charmed mesons (D, and D) with light mesons {n and K). In reference [6], 
Dn and DK scattering were studied using the scalar form factors in semileptonic pseudoscalar-to- 
pseudoscalar decays. In our work we use Liischer's formula [7] to extract scattering information 
from the energy shift of two interacting hadrons relative to the total energy of the two individual 
hadrons. 



2. Fermion Actions 

We use Fermilab formulation [8] for the charm quark, domain-wall fermions for the light 
quarks and staggered sea quarks. In the following, we will specify the Fermilab formulation as 
well as the tuning of the parameters in this formulation. The action is: 

S = Sq + Sb + Se, 

b b 
Sq = £^(x)[mo + (7oVo--Ao) + v£(T4V,---A,)]?W, 

b 

Sb = --CBY,^{x)(^(^ijFij)q{x), 

^ .V i<j 

b 

Se = --c£^^(.x)(^ao,Fo,-)^(.x), 

^ X i 

where b is the lattice spacing, Vq and V; are first-order lattice derivatives in time direction and 
space directions, Aq and A,- are second-order lattice derivatives, F^v is the field tensor defined 
in reference [8]. In this action, the space-time exchange symmetry is not imposed. Sq is just 
the standard Wilson fermion action except that the coefficient in front of the space term and the 
coefficient in front of the temporal term are different. Sb and Se are spatial and temporal clover 
terms, also with different coefficients. 

There are four pai^ameters to tune: the charm quark mass m^, the anisotropy v, and the two 
clover coefficients cg and ce- We use the spin-average mass of and Tj^ to tune the charm quark 
mass. The value of v is tuned to restore the dispersion relation. As for cb and ce, the tree level 
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tadpole improvement estimate is cb = ce - 
parameters [9]: 

V 

Cb = 

Here cb and ce both depend on v. We use 

3. Numerical Ensembles 

We employ the gauge configurations generated by the MILC Collaboration [10]. We use the 
20^ X 64 lattices generated at four values of Ught-quark masses. The lattice spacing b = 0.12406 
fm. The details of the ensembles are listed below: 



Ensemble 


bnii 




bmj^^'^ 


bnif"'^ 


number of props 


2064f21b676m007m050 


0.007 


0.050 


0.0081 


0.081 


450 


2064f21b676m010m050 


0.010 


0.050 


0.0138 


0.081 


650 


2064f21b679m020m050 


0.020 


0.050 


0.0313 


0.081 


550 


2064f21b781m030m050 


0.030 


0.050 


0.0478 


0.081 


380 



The subscript / denotes light quark, and s denotes the strange quai^k. The superscript dwf denotes 
domain-wall fermion. 

4. Heavy- Quark Action Test 

In the heavy-quark action, the dispersion relation f"^ = + c^p^ needs to be restored by tun- 
ing the parameter v to get the value of to agree with the theoretical value 1 . To do that, we calcu- 
late the single-particle energy of rjc, Dg and D at the six lowest momenta: ^ (0, 0, 0) , ^ ( 1 , 0, 0) , 
^ ( 1 , 1 , 0) , ^ ( 1 , 1 , 1 ) , ^ (2, 0, 0) , ^ (2, 1 , 0) . We tune the dispersion relation of Tj,. and get the dis- 
persion relations of and D to be restored as well. The following table lists the values of 
we get from the numerical simulation: 







J/'¥ 


D 




c2 


0.989(0.005) 


0.965(0.009) 


1.012(0.017) 


1.006(0.009) 



By tuning the charm-quark mass, we get the spin-average mass of rjc and7/*P to be 3056.54(1.15)MeV, 
which agree well with the experimental value. We also calculate the hyperfine splitting, which we 
find to be 99.1(l.l)Mev. 

In this section, all numerical results are obtained on the ensemble 2064f21b676m007m050 
which has the lightest light-quark mass. We keep the charm quark mass and anisotropy parameter 
fixed for all ensembles. 

5. Baryon Spectrum 

We calculate the masses of singly charmed baiyons and doubly charmed baryons at four differ- 
ent light quark masses. We employ a simple linear- relation rtiharyon = c\ +C2m^//| to extrapolate 
the baryon masses to the physical point. The values of mj^/ fji are taken from reference [11]. Fig. 
1 shows the baryon mass spectrum. We present the singly charmed baryon mass splittings in Fig. 
2. Our results ai^e compai^able with similar work done using a staggered light-quark a ction [12]. 



1/mq. Chen suggested a better way to evaluate the two 
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Chen's evaluation in oui" work. 



CE = ^{l + v)^. (2.1) 
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Figure 1: Charmed baryon masses as functions of m^/ f^. The masses are divided by /o = 130.7MeV to 
make them dimensionless. The four blue points in each panel denote the baryon masses at four different 
light-quark masses. They are extrapolated to the physical values, which are denoted by the red points. 
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Figure 2: Mass splittings of singly charmed baryons. The red lines represent the experimental values. The 
four bars in different colors represent the numerical values of the mass difference between A^- and other 
singly charmed baryons. The heights of the bars indicate the statistical errors. 



6. Charmed Hadron Interactions 

Liischer has shown that the scattering phase shift is related to the energy shift {AE) of two 
interacting hadrons relative to the total energy of the two individual hadrons [7]. The total energy 
of two interacting hadrons (hi and /j2) is obtained from the four-point correlation function: 



G 



t\ 



(6.1) 
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To extract the energy shift /\E, we define a ratio ^^{t): 

R''''^-{t) = ^, — , V X ^exp(-A£-0, (6.2) 

where G''' (f,0) and G^^{t,0) are two-point functions. AE can be obtained by fitting R'^^^^^{t) to a 
single exponential. 

The magnitude of center-of-mass momentum p is related to /\E by 



/\E = ^ ^m\^^ ^ ^ m\^_ - ruhi - ruh^ . (6.3) 
The phase shift is obtained from the following relation: 

pc„.5(p) = l-s((|)=). (64) 

where the S function is defined as 

jl<A 1 

S(x)= £ 47rA. (6.5) 

The sum is over all three-vectors of integers j such that | j | < A, and the limit A ^ oo is implicit. If 
the interaction range is smaller than half of the lattice size, the 5-wave phase shift will be 

pcot5(p) = ^ + ^(p2), (6.6) 

where a is the scattering length. 

By measuring the energy shift, the momentum p can be obtained by equation (6.3). Then we 
calculate the right-hand side of equation (6.4). From equation (6.4) and equation (6.6), we can see 
that the scattering length is just the inverse of the right-hand side of equation (6.4). 

We study the scattering of the charmonia (t]c, //*!') with light hadrons {n, p, N). The scattering 
of charmed mesons {Ds,D) with n and K are also studied. We choose the isospin-| channel for 
D — n scattering and isospin-1 channel for D — ^ scattering to avoid disconnected diagrams. For 
the same reason, we use D+ and for Dg — K scattering. The interactions of the heavy hadron(/j) 
with the pion has a special feature due to the Nambu-Goldstone nature of the pion. The 5-wave 
scattering length is given by [13] 

^h-n ^ ^ i) _;^(;,, + 1) -2] + i?{ml), (6.7) 

where I is the total isospin of the 7r — /i system. The first term in equation (6.7) vanishes for — tt, 
Tjc — 71 and D.V — 7i channels but not for D — tt channel. So the leading term of D — tt scattering 
length is proportional to niji/ fji while the leading terms of — n, rjc — 71 and — 7i scattering 
lengths are proportional to m^/f^ ■ We measure the scattering lengths at four different light-quark 
masses. We use the relation a = ci+ cjmji/fji to extrapolate D — 7i scattering length to the physical 
point. For other channels, we use the relation a = ci +C2?m^//| to do the extrapolations because 
the leading term in equation (6.7) vanishes. Fig. 3 and Fig. 4 show the scattering lengths of all 
these channels. 
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Figure 3: Charmonium to pion, rho and nucleon scattering lengths as functions of rn^/ f^. In each panel, the 
four blue points denote the scattering length measured at four different light-quark masses. The red points 
denote the extrapolated values. 




Q 

CO 




Figure 4: Charmed meson scattering lengths as functions of m^/ f'^{m-,i/ fn). In each panel, the four blue 
points denote the scattering length measured at four different light-quark masses. The red points denote the 
extrapolated values. 



The following tables list the extrapolated scattering lengths of all channels: 



Chaimel 


Scattering lengths(fm) 


Channel 


Scattering lengths(fm) 


77c -71 


0.00(1) 


7/*P - p spin 1 


0.01(6) 




0.03(5) 


J/^-p spin 2 


0.01(6) 


ric-N 


0.18(9) 


D-K 


-0.16(4) 


J/^'-n 


-0.01(1) 


D-K 


-0.23(4) 


J/^-N spin 1/2 


-0.05(77) 




0.00(1) 


J/^'-N spin 3/2 


0.24(35) 


D,-K 


-0.31(2) 


J/^-p spinO 


0.00(6) 
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7. Conclusion 

We calculated the scattering lengths of charmed mesons and charmonia with light hadrons. 
For the channels of charmonia with light hadrons and — 71 channel, we found weak interactions. 
The scattering lengths are zero or close to zero. For \he D — n^D — K and Dg — K chaimels, we 
found relatively strong repulsive interactions. In the future, we will improve statistics to get more 
reliable results, and perform extrapolations using chiral perturbation theory forms. 
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